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We update our analysis of p-ui mixing effects in the pion form factor to incorporate 
the recently published BaBar e + e~ — > ir + ir~ cross-sections. The implications for r- 
decay-based Standard Model estimates of the leading order hadronic contribution, 
l a dhad> t° the anomalous magnetic moment of the muon, and for the extraction of 
the off-diagonal vector meson self-energy matrix element, n 



m 



are discussed. 



PACS numbers: 13.66.Bc,13.75.Cs,14.60.Ef,13.40.Em 



In the following we update the analysis performed in Ref. [1| of the isospin-breaking 
(IB) p — bj mixing correction required in order to use r-decay-based data instead of 
electroproduction data in the evaluation of [a^]^' the leading order hadronic vacuum 
polarization contribution to the anomalous magnetic moment of the muon, = (g — 
2)^/2. This update focusses on the BaBar electroproduction data 0, |3[ since it was not 
released until shortly after publication of the previous analysis. 

As is well known, several recent measurements of the e + e~ — > 7r + 7r~ cross-section 
|9[ together yield estimates of [a^]^ which are consistent with one another, but lead to 
Standard Model (SM) predictions for a M deviating from the BNL E821 experimental 
result To|| by ~ 3.2 — 3.6a |11 - 18| . In contrast, evaluating [a^jh^d usm g T decay data in 



place of isovector electroproduction data 
from experiment by only ~ 1.9 — 2.4<r 11 



19 



12 



24| yields a SM prediction for a M differing 



15Hl8|. Use of the r decay data requires 



that a number of small IB corrections to the CVC relation be taken into account. These 



corrections have been extensively studied in Refs. [16|, l25H29| and are believed to be well 
understood. We denote these corrections collectively by [^a^]^ an d focus in what follows 
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on the particular contribution arising from p — u mixing, [Sa p ]^ d . mix , which is defined 
explicitly in Ref. 



An important observation made in Ref. [30| was that the generic structure of the 
p — bj interference contribution to F^{s) introduces strong fit-parameter-sensitive can- 
cellations, and hence significant model dependence, into the integral corresponding to 
[^ a fj]had mix- O ur analysis thus employs a range of models for F^(s),all having some basis 
in phenomenology These are the Kuhn-Santamaria (KS) model |32Li_the Hidden Local 
Symmetry (HLS) model (H, 34 L the Gounaris-Sakurai (GS) model 35|, and a modified 
version of the GP/CEN model 26 L 36|. (Detailed descriptions of the models can be found 
in section II of Ref. Q.) Refs. IT 13 01 ] show that it is necessary to consider such a range of 
models if one wishes to properly assess the model dependence of [Sa^]^ d . mix , and, from 
this, the uncertainty in the nn contribution to the r decay-based estimates of [aJhS- 

Shortly after the publication of Ref. [l[ the BaBar collaboration released the data 
corresponding to its measurement of the e + e~ — > 7r + 7r~(7) cross-section, using the initial- 
state radiation method, from threshold to 3 GeV 0, Compared to the electroproduc- 
tion data sets described and used in Ref. Q (CMD-2 SND and KLOE 0,0) 

the BaBar data offers considerably increased statistics, including 15 data points in the 
interference region (770-800 MeV), as well as generally lower statistical and systematic 
errors. The BaBar data distinguishes itself from its predecessors, however, in that the 
value of a p computed using it as the source of the im contribution to [fl/Jhad more closely 
corresponds to the experimental and r decay based values, deviating from the experi- 
mental value by only 2.4<r (37[. 

As before, we perform fits to the BaBar data set using the models indicated above. 
Although the BaBar data extends up to 3 GeV, only the low-energy part of this data is 
relevant to analyzing p — u mixing. We, therefore, limit our analysis to the maximum 
e + e~ center-of-mass energy of 970 MeV employed in our previous analysis. The results 
quoted below for [Sa p }^ d . mix are insensitive to modest changes in this choice of endpoint. 
All results correspond to the bare form factor (i.e. with the effects of vacuum polarization 
removed). Details of the fit procedure, including all input values, are unchanged from 
Ref. Fit results for each model are shown in Table [H The fit parameters are the p mass 
and width, m p and T p , the complex coefficient of the u> contribution, 5, the coefficient of 
the p' term, 0, and the HLS model parameter, cihls- A blank entry indicates that a fit 
parameter is inapplicable to that particular model. For the GP/CEN + and GP/CEN ++ 
models, the effective value of T p is shown in brackets to highlight that it is in fact 5T p , an 
offset from the nominal chiral effective theory p width, which is the actual fit parameter. 

Comparing the results of Tabled with those of Tables I-IV in Ref. we see that the 
BaBar data yields a p width larger by 1-7 MeV (depending on the specific data set and 
model) and a reduced p-uj mixing phase. Reasonable x 2 /dof results are obtained despite 
the reduced scale of statistical errors in the BaBar data relative to the other data sets. 

The values obtained for [Sa p ]^ d . mix using the BaBar data and for each of the models 
considered are shown in Table [Hi along with the values from the data sets used in Ref. [l[. 
The latter are included for ease of comparison. The BaBar data yields somewhat larger 
central values, along with reduced errors. However, as before [J, 3(|, the variation in 
the values of [<5cv]£^. m j X across the various models is greater than the experimental 
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TABLE I: Results of fits to the BaBar 2009 data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


m p (MeV) 


i 1 i r\ on 

772.11±0.30 


773.48±0.29 


'-j'i A on i n on 

774.29±0.30 


T^7f o'i i n nn 

775.87±0.29 


775.87±0.29 


T p (MeV) 


147.56±0.54 


149.68±0.57 


149.87± 0.57 


(148.66) 


(148.65) 


5T p (MeV) 








1.33±0.46 


1.32± 0.46 


\5\ (10- 3 ) 


1.89±0.03 


1.99±0.03 


1.96±0.03 


2.30±0. 03 


1.98± 0.03 


Arg(<5) (deg) 


9.4 ± 1.2 


10.3 ± 1.1 


10.3 ± 1.1 


10.8 ± 1.1 


10.8 ± 1.1 


P 


-0.152±0.002 




-0.088±0.002 






O-HLS 




2.3989±0.007 








X 2 /dof 


392/250 


320/250 


322/250 


413/251 


414/251 



TABLE II: [ba p ]f° d . mix x 10 10 for the models discussed in the text and the CMD-2, SND, and 
KLOE e + e~ — > it + tt~ cross-sections. 



Experiment 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


CMD-2(94) 


3.8 ±0.6 


4.0 ±0.6 


2.0 ±0.5 


2.0 ±0.5 


1.8 ±0.4 


CMD-2(98) 


4.0 ±0.6 


4.6 ±0.6 


2.5 ±0.5 


2.2 ±0.4 


2.1 ±0.4 


SND 


4.2 ±0.4 


4.3 ±0.4 


2.2 ±0.3 


1.9 ±0.3 


1.7 ±0.3 


KLOE(02) 


(2.2 ±0.6) 


4.2 ±0.7 


2.2 ±0.6 


(0.5 ±0.8) 


(0.3 ±0.8) 


BaBar(09) 


5.0 ±0.2 


5.0 ±0.2 


2.9 ±0.2 


2.6 ±0.2 


2.4 ±0.2 



uncertainty produced by any single model. In arriving at a final assessment of our results 
for [8 a n] hod-mix' we nave adopted the view that, since all the models considered have 
a reasonable basis in phenomenology, all results corresponding to a given data set and 
given model which produce an acceptable quality fit are to be included in the assessment. 
(Those entries in brackets in Table ILT1 correspond to poor quality fits and are not included 
in our final result.) We thus first perform a weighted average over all experiments for 
each separate model, and then take the average (half the difference) of the maximum and 
minimum values allowed by the resulting error intervals for the different models to define 
our central values (model-dependence-induced uncertainties). The updated combined 
assessment, now including the BaBar results, is 

[H]L' ;mra = (3-5 ± l.5 model ± 0.2 data ) x 1(T 10 . (1) 

The central value has increased by 0.4 x 10 -10 and the data error has decreased by 
0.1 x 10~ 10 compared to the value reported in 

The value shown in Table ILT1 obtained using BaBar data and the GS model is compat- 
ible with the GS result reported in Ref. j37|- The KS model result, however, is not, the 
KS and GS results for [^a^jf^ . differing significantly in Table ILT1 but being the same 



in Ref. [37|. The source of this apparent discrepancy is that two distinct 'KS' models 



have in fact been employed: the one we denoted KS above, and the alternate version 
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used in Ref. [371, which we call KS'. As discussed in Ref. Ill these two models differ in 



the s-dependence assumed for the p-u mixing contribution to F n (s). We have confirmed 

iLO 

W hadimix 



that the alternate, KS', form indeed yields results for [Sa^] had , mix compatible with those 



of the GS model. In fact, it turns out that the presence or absence of the extra s/mf, 
factor (which is what distinguishes the KS and KS' model forms) is also the key feature 
distinguishing those models which yield 'high' values of [&fy]f a< f. mia . (KS, HLS) from those 
which yield 'low' values (GS, GP/CEN ++ ). The data, in the narrow range of s over which 
p — bj interference is significant, is incapable of distinguishing between these differing s- 
dependences. While such differences have only a very small impact on the values of the 
model fit parameters, the presence or absence of the factor of s/m^ strongly affects the 
very close cancellation occurring in the weighted integral for [5a^]f d ■ ■ Since there is, 
at present, no compelling theoretical argument favouring one choice of s-dependence over 
the other in the interference region, we adopt the view that the unknown s-dependence 
of the mixing term must be treated as an additional source of uncertainty for [5a M ], , . . 
This uncertainty significantly increases the total error on [a ^ad-mix- 

As explained in Refs. [l|, 38j, analysis of the electroproduction data in the interference 
region also allows one to extract the off-diagonal p — u element of the vector meson 
self-energy matrix, U pLJ (q 2 ), and the isospin-breaking coupling ratio G = g^i^j 'sy^, 
with g^i-K-K and g p i n7r the isospin-pure tttc couplings of the p and uj mesons. Ii pul {q 2 ) is 
of interest, for example, for meson-exchange models of IB in the NN interaction. The 
procedure for performing this determination has been described in detail in Refs. fll. |38|. 

The separation of mixing and direct u> — > ttti contributions depends on the model 
used for the broad p contribution to F n (s). We report in Table Hill the results for 
(the Orsay phase), G, and T = lJ puj (m 2 p )/rh p T p (with Il puJ the real part of U puj and rh p 
the real part of the complex p pole position), obtained from the BaBar data set for the 
various models used. The one-sigma contours for G and T are shown in Fig. [TJ The 
corresponding contours for the CMD-2(98) and SND data sets are shown for comparison 
in Fig. [2j Readers are directed to Ref. |l( for full details. 



TABLE III: Orsay phase and separated mixing and direct ojtttt coupling parameters for the 
BaBar(09) data. 



Parameter 


KS 


HLS 


GS 


GP/CEN+ 


GP/CEN++ 


<P (deg) 


108 ±1 


108 ±1 


107 ±1 


107 ±1 


107 ±1 


G 


0.028 ± 0.013 


0.035 ± 0.013 


0.036 ±0.013 


0.039 ±0.015 


0.040 ± 0.013 


f 


-0.037 ± 0.002 


-0.038 ± 0.002 


-0.038 ± 0.002 


-0.0381 ± 0.0009 


-0.038 ± 0.001 



The results of Table HTT1 should be compared to those in Tables VI to IX of Ref. [lj. It is 
immediately apparent that the lower statistical uncertainty of the BaBar data translates 
into much greater precision in the extracted value of T. Two further significant differences 
concern the central values of <fi and G, which are both lower for the BaBar data compared 
to the other data sets. The BaBar data also significantly improves the significance of 
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FIG. 1: BaBar(09) G and T one-sigma re- FIG. 2: CMD-2(98) and SND G and T one- 
gions. sigma regions. 

the evidence for G ^ 0. In Ref. [l[ we presented combined averages both including the 
KLOE data and excluding it. The high precision BaBar data now so dominates the 
combined averages that there is little distinction between the results obtained including 
or excluding KLOE; we thus present only the former in Table HVl below. 



TABLE IV: Combined averages including/excluding BaBar data. 





This work 


Ref. [1] 


KLOE Included 




109.0° ±1.9^ ±0.8^ 


1 1 QO 1 A O 4-9° 


f 


-0.041 ± 0.003 mO(fei ± 0.001 data 


-0.044 ± 0.006 modei ± 0.002 data 


G 


0.054 ± 0.014 modei ± 0M0 data 


G = 0.080 ± 0.026 model ± 0M5 data 



Note that the lower model dependence shown in the first column of Table IIVI reflects 
the dominance of the high-precision BaBar data over the other data sets in the averages, 
rather than any improved model consistency. The combined average for the complex- 
valued off-diagonal part of the physical p — u self-energy matrix, U pLU (m 2 p ), now including 
the BaBar data, is 

iWmJ) = (-4620 ± 220 modcl ± 170 data ) + (-6100 ± 1800 modcl ± 1110 data )z MeV 2 . (2) 

In conclusion, we have updated the determination of [(Jo.^] . . . and the separation of 
p—uo interference in the e + e~ — > 7t + tt~ cross-sections into direct and mixing induced terms 
using the recently released BaBar ISR data. The main results are given in Eqs. [TJ and HI 
and in Table [TTT1 We conclude that, while not at present dominant, the model-dependence 
of [$ a ij]h d- m ix given in Eq. [TJ will eventually represent a fundamental limitation on the 
use of t data in the evaluation of a M . 
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